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2 ■ Abstract 

We analyze the proton and deuteron structure functions at large x using a re- 
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CO , cently introduced scaling variable x. This variable includes power corrections 

CN 

to x-scaling, and thus allows us to reach the Bjorken limit at moderate Q 2 . 
\ Using available data we extract the ratio F^ip^/F^ix) for x < 0.85. Contrary 
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to earlier expectations this ratio tends to the value ~2/3 for x — > 1, which 
corresponds to the quark model prediction for equal distributions of valence 



According to the quark-parton model, the structure functions of hadrons in the Bjorken 
limit (Q 2 = q 2 — v 2 — > oo and x = Q 2 /2iWV=const) are directly related to the parton 
distributions qi(x). For instance 

F 2 (x,Q 2 ) -> F 2 {x) = J £e 2 i xq i {x), (1) 

i 

where the sum is over the partons, whose charges are e,. In the region x — > 1, the contribution 
of sea quarks can be neglected. Then assuming the same distribution of the valence quarks, 
one easily finds from Eq. (|XJ) that the neutron-to-proton ratio F^x) / F 2 {x) approaches 2/3 
for x — > 1. If, however, the quark distributions are different, one can establish only upper 
and lower limits for this ratio, 1/4 < F£ / ' F% < 4, which follow from isospin invariance 



IJ. The existing data show considerable Q 2 -dependence of the structure functions which is 
attributed mainly to the QCD logarithmic corrections to Bjorken scaling. However, at high- 
x the scaling violations are dominated by power corrections oc 1/Q 2 (higher twist and target 
mass effects), which are difficult to evaluate. Therefore, in order to check the parton model 
predictions, one needs to obtain the structure functions at high Q 2 , where these corrections 
are small. At present, high Q 2 structure functions (Q 2 ~ 250 (GeV/c) 2 ) extracted from 
BCDMS [§ and NMC [|] data are available only for x < 0.7. The ratio F£{x)/F%(x) 
obtained from an analysis of these data f||| shows steady decrease with x. Thus, it is 
usually assumed that this ratio would reach its lower bound, F^jF^ —>■ 1/4, for x — > 1. 
This corresponds to d(x)/u(x) — > for x — > 1, where d(x) and u(x) are the distribution 
functions for up and down quarks in the proton. To check this assumption one needs data 
for the structure functions at x <^ 0.7. The latter are now available only at moderate values 
of momentum transfer, Q 2 < 30 (GeV/c) 2 f6|-f|, and exhibit very strong Q 2 -dependence. 

A part of the Q 2 -dependence of the structure functions, generated by the target mass 
corrections, is usually accounted for by using the Nachtmann scaling variable pH|jn[ ] 
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instead of the Bjorken variable x. The question is whether the target mass effect is respon- 
sible for a major part of the scaling violation at high x. If so, the replacement of x by £ 
in structure functions would allow us to reach the scaling limit already at moderate Q 2 . 
However, the analysis of recent high x SLAC data in terms of the Nachtmann variable 
still reveals strong Q 2 -dependence of the structure functions. 

Besides the target mass effects, one can expect important nonperturbative effects from 
the confining interaction of the partons in the final state. Indeed, the partons are never 
free, so that the system possesses a discrete spectrum in the final state. Although in the 
Bjorken limit the struck quark can be considered a free particle, the discreteness of the 
spectrum manifests itself in 1/Q 2 corrections to asymptotic structure functions JT2]ir3] ]. One 



can expect that these corrections are significant in particular at high x, where lower-lying 



excitations should play an important role. 

We have found jnj that the target mass and confining interaction effects in the final 
state can be effectively accounted for by taken the struck quark to be off-shell, with the 
same virtual mass before and after the virtual photon absorption. As a result, the Bjorken 
scaling variable x is replaced by a new scaling variable x = x(x, Q 2 ), which is the light-cone 
fraction of the off-shell struck quark. Explicitly, 



x + y/1 + AM 2 x 2 /Q 2 - - x) 2 + 4m 2 x 2 /Q 2 

1 + JT+ 4M 2 x 2 /Q 2 ' ^ 



where M is the target mass and m s is the invariant mass of spectator partons (quarks and 
gluons). For Q 2 —> oo or for x — > the variable x coincides with the Nachtmann variable £, 
Eq. (Q). However, at finite Q 2 these variables are quite different. 

It follows from Eq. (§) that x depends on the invariant spectator mass, m s . The latter 



can be considered a function of the external parameters only In the limit x — > 1 

(elastic scattering) no gluons are emitted, and thus m s — > m , the mass of a two-quark 
system (diquark). When x < 1, the spectator mass m s increases due to gluon emission. For 
x close to 1, ra s can be approximated as 



m 



m 2 + C(l-x), (4) 



where the coefficient C ~ (GeV) 2 . In the following we regard it as a phenomenological 
parameter, determined from the data. 

Consider the proton and deuteron structure functions (per nucleon) ,d (x, Q 2 ) at large 
x. These are shown in Fig. la,b as functions of x for Q 2 =230 (GeV/c) 2 , which is the 
maximal value of Q 2 in the BCDMS measurements 0. (Notice that for such high values 
of Q 2 the variable x is close to x, Eq. ([$])). The data are taken from 0, where the solid 
lines correspond to a 15 parameter fit to BCDMS and NMC data 0. For smaller values of 
Q 2 , the violations of Bjorken scaling are very significant, Fig. lc,d. By assuming that the 
major part of x-scaling violations at high x are correctly accounted for by the variable x, 
the Q 2 -dependence of the structure functions in this region is given by 



FZ> d (x,Q 2 ) = F? d (x(x,Q 2 )). (5) 

Let us compare this result with two BCDMS and SLAC data bins for x=0.65, 0.75 Fig. 
lc,d, which are the largest values of x available in the BCDMS experiment. We find that 
these data are perfectly reproduced (the dashed lines in Fig. lc,d), by taking the spectator 
mass m 2 =0.75 (GeV) 2 for x=0.75 and m 2 =1.05 (GeV) 2 for a;=0.65. It is quite remarkable 
that the same values of m s are obtained for proton and deuteron targets, although the 
corresponding structure functions are rather different. Using Eq. one finds that these 
values of m s determine the parameters mo and C, namely C = 3 (GeV) 2 and mo = 0. The 
latter implies that the spectator quarks are massless and collinear. 

Using these values of mo and C for definition of m s in Eq. @ we can study the structure 
functions for x > 0.75. Consider first the data for the proton structure function from 
the SLAC experiments in the region x > 0.6 for 5 ^ Q 2 ^ 30 (GeV/c) 2 , plotted 

as a function of £ and x respectively, Fig. 2a,b. The data points close to the region of 
resonances were excluded by a requirement on the invariant mass of the final state, namely 
(M + v) 2 -q 2 > (M + A) 2 , where A =300 MeV. The solid line and the three data points in 
Fig. 2a,b correspond to the asymptotic structure function at Q 2 = 230 (GeV/c) 2 , the same as 
in Fig. la. The analysis in terms of the Nachtmann scaling variable £, Fig. 2a, shows poor 
scaling. Moreover, the data points are far off the asymptotic structure function (the solid 
line). In contrast, the same data plotted as a function of x, Fig. 2b, show excellent scaling. 
Also the data points completely coincide with the structure function at Q 2 = 230 (GeV/c) 2 , 
available for x < 0.75. This agreement provides strong evidence that the x-scaling is not 
accidental. We therefore propose that the data points in Fig. 2b represent a measurement 
of the asymptotic structure function for x > 0.75 as well. 

Next, consider the deuteron structure function from the SLAC data 0-|8[], Fig. 3a,b. 
As in the previous case, we exclude the region of resonances by taking the invariant mass 
in the final state greater then M + A, where A=300 MeV. In addition, in order to avoid 
complications from binding and Fermi motion effects, we exclude from our analysis the data 
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points with x > 0.9. (This restriction is relevant only for the data ||). Indeed, recent 
calculations of Melnitchouk et al. |T3] show that the ratio IF^jiF^ + F£) is about 1.13 for 
x = 0.9 and Q 2 =5 (GeV/c) 2 , and it rapidly increases for x > 0.9. However, for x < 0.85, 



this ratio is within 5% of unity Fig. 3a shows the deuteron data plotted as a function of 
the Nachtmann variable £. The solid line and three data points show the structure function 
at Q 2 = 230 (GeV/c) 2 , the same as in Fig. lb. One finds that the data display no scaling 
and they are far from the solid line. The same data as a function of x are shown in Fig. 3b. 
As in the proton case the data show very good scaling and do coincide with the structure 
function at Q 2 = 230 (GeV/c) 2 . Unfortunately, there are no deuteron data in the high-x 
region for large Q 2 , as for instance the proton data 0. As a result, the available deuteron 
data allow us to determine the asymptotic structure function only up to x = 0.85. 

Now with the asymptotic structure functions F 2 (x) = F 2 (x) = F 2 (x,Q 2 — > oo) found 
above, we can obtain the ratio F 2 (x) / F 2 (x) . For this purpose we parametrize the asymptotic 
structure functions as F 2 d {x) = exp(— J2t=o a-iX 1 ) and determine the parameters a; from the 
best fit to the data in Figs. 2b, 3b. The resulting F%jF 2 ratio is shown in Fig. 4 by the solid 
line. The dotted lines are the error bars on the fit, which combine statistical and systematic 
uncertainties. The dashed line corresponds to F 2 /F 2 = 2/3. For a comparison, we show 
by the dot-dashed line a polynomial extrapolation of this ratio to large x, obtained from 
BCDMS and NMC data by assuming that F 2 n /F 2 p -> 1/4 for x -> 1 §. 

Our results shown in Fig. 4 demonstrate that contrary to earlier expectations, the ratio 
F 2 /F 2 does not approach its lower bound, but increases up to approximately 2/3. The 
latter is the quark model prediction, assuming identical distributions for each of the valence 
quarks. The accuracy of our results will be checked in future experiments, which will provide 
high Q 2 data for the structure functions at large x. 
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FIGURES 

FIG. 1. (a,b) F%> d (x) = F%> d (x,Q 2 ) at Q 2 =230 (GeV/c) 2 . The data are from BCDMS mea- 
surements [g]. The solid line is the 15 parameter fit [|| to BCDMS, NMC and SLAC data. (c,d) 
Proton and deuteron structure functions at constant x. The dashed lines show Q 2 -dependence 
of the structure functions given by Eq. (||). The data are from BCDMS || and SLAC HQ] 
experiments. The error bars show combined statistical and systematic errors. 

FIG. 2. (a) SLAC data for the proton structure function for 5 < Q 2 < 30 (GeV/c) 2 , plotted as 
a function of the Nachtmann variable £, Eq. (||). The data with largest value of £ are taken from 
recent measurements fj|. Three high-statistics data sets for Q 2 =5.9, 7.9, and 9.8 (GeV/c) 2 , taken 
from HH, are marked by "+", "x", and respectively. The other data points are from Hffj. 
The error bars show combined statistical and systematic uncertainties. Three data points marked 
by "o" and the solid curve are the same as in Fig. la, and show the structure function at Q 2 =230 
(GeV/c) 2 . (b). The same data plotted as a function of the scaling variable x. The data points for 
5 < Q 2 < 30 (GeV/c) 2 coincide well with the structure function at Q 2 =230 (GeV/c) 2 . 

FIG. 3. (a) SLAC data for deuteron structure function for 5 < Q 2 < 30 (GeV/c) 2 plotted as a 
function of the Nachtmann variable £, Eq. (g). Four high-statistics data sets for Q 2 =3.9, 5.9, 7.9, 
and 9.8 (GeV/c) 2 , taken from [§, are marked by "*", "+", "x" and "#" respectively. The other 
data points are from The error bars show combined statistical and systematic uncertainties. 

Three data points marked by "o" and the solid curve are the same as in Fig. la, and show the 
structure function at Q 2 =230 (GeV/c) 2 . (b). The same data plotted as a function of the scaling 
variable x. 

FIG. 4. Neutron-to-proton structure function ratio at large x. The solid line is the result of 
our analysis. The dotted lines show combined statistical and systematic errors. The dashed line 
is the quark model prediction for x — > 1 for equal distributions of valence quarks. The dot-dashed 
line shows the expected behavior of this ratio from polynomial extrapolation of BCDMS and NMC 
data |j . 
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